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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

meters 

gallons  (US  liquid) 

3.785412 

cubic  decimeters 
(litres ) 

inches 

2.54 

centimeters 

kips  (force) 

4.448222 

kilonewtons 

kips  (force)  per 

6.894757 

megapascals 

square  i'nch 

megatons  (nuclear 

4.184 

petajoule3 

equivalent  of  TNT) 

pounds  (force)  per 

6.894757 

kilopascals 

square  inch 

pounds  (mass) 

0.4535924 

kilograms 

pounds  (mass)  per  cubic 
foot 

16.01846 

kilograms  per  cubic 
meter 

DEVELOPMENT  OF  A  THEORETICALLY  SOUND  CAP  MODEL 
FOR  FITTING  ISST-TYPE  MATERIAL  BEHAVIOR 


INTRODUCTION 

Incremental  elastic-plastic  constitutive  models  which  include  both  an 
ultimate  failure  envelope  and  a  work-hardening  yield  surface  —  or  cap  —  are 
commonly  used  in  ground  shock  calculations  to  simulate  geologic  material 
behavior  (Reference  1).  The  elastic  3hear  modulus  in  the  cap  model  is  usually 
formulated  in  terms  of  the  second  invariant  of  the  stress  deviation  tensor 

p 

JJ  and  the  plastic  volumetric  strain  ekk  .  Such  formulations,  however, 
cannot  replicate  the  highly  nonlinear  bowl-shaped  unloading  stress  paths  that 
typify  ISST-type  material  behavior  (Reference  2  and  Figure  1). 

In  1 98M ,  a  cap  model  was  developed  in  which  the  shear  modulus  was  formu- 

P 

lated  not  only  as  a  function  of  JJ  and  ekk  ,  but  also  as  function  of  the 

third  invariant  of  the  stress  deviation  tensor  JJ  (Reference  3).  Although 

this  model  fit  the  recommended  ISST  material  properties  quite  well,  it  did  not 

satisfy  all  of  the  vigorous  theoretical  constraints  outlined  by  Dr.  Ivan  S. 

Sandler  of  Weidlinger  Associates  as  a  prerequisite  for  use  of  a  constitutive 

model  in  explosive-induced  ground  shock  and  soil/structure  interaction  codes 

(Reference  4).  In  order  to  satisfy  these  theoretical  requirements  and  3till 

satisfactorily  replicate  ISST-type  material  behavior,  a  new  cap  model  has  been 

p 

formulated  in  which  the  shear  modulus  is  expressed  as  a  function  of  ekk  and 
the  second  invariant  of  the  strain  deviation  tensor  IJ  .  In  addition  to  its 
theoretical  soundness,  the  model  has  the  advantage  for  finite  element  calcula¬ 
tions  of  not  requiring  any  additional  storage  since  IJ  can  be  calculated 
directly  from  displacements.  Strain-dependent  models  also  have  an  advantage 
over  stress-dependent  models  in  that  they  can  be  numerically  incorporated  into 
computer  codes  so  as  to  be  insensitive  to  increment  size. 

A  general  description  of  the  strain-dependent  cap  model  and  the  associ¬ 
ated  theoretical  proofs  follows. 

GENERAL  DESCRIPTION  OF  THE  NEW  MODEL 

The  basic  premise  of  elastic-plastic  constitutive  models  is  the  assump¬ 
tion  that  certain  materials  are  capable  of  undergoing  small  plastic  (perma¬ 
nent)  as  well  as  elastic  (recoverable)  strains  at  each  loading  increment. 


Mathematically,  the  total  strain  increment  is  assumed  to  be  the  sum  of  the 
elastic  and  plastic  strain  increments;  i.e.. 


deu  ■  “'fj  *  defj 


where 

de^  =*  components  of  the  total  strain  increment  tensor 

g 

dEjj  -  components  of  the  elastic  strain  increment  tensor 

p 

de^j  *»  components  of  the  plastic  strain  increment  tensor 

Within  the  elastic  range,  the  behavior  of  the  material  can  be  described 
by  an  elastic  constitutive  relation  of  the  type 


de.  .  =  C.  ,  (o  )  do,  , 
lj  ljkl  mn  kl 


where 

cljkl  "  material  response  function 

dokl  -  components  of  stress  increment  tensor 

The  behavior  of  the  material  in  the  plastic  range  can  be  described  within  the 
framework  of  the  generalized  incremental  theory  of  plasticity.  The  mathemati 
cal  basis  of  the  theory  was  established  by  Drucker  (Reference  5),  who  intro¬ 
duced  the  concept  of  material  stability;  this  concept  has  the  following 
implications: 

1.  The  yield  surface  (loading  function)  should  be  convex  in  stress 
space. 

2.  The  yield  surface  and  the  plastic  potential  should  coincide  (which 
results  in  an  "associated"  flow  rule). 

3.  Work  "softening"  should  not  occur. 

These  three  conditions  can  be  summarized  mathematically  by  the  following 
inequality 


Elastic  Strain  Increment  Tensor 

1  > 

For  isotropic  elastic  materials,  the  strain  increment  tensor  (Equation  2) 
takes  the  following  form 


E  dJ1  1 

deij  “  9K  6iJ  +  2G  dSij 
or 

d»u  ■  K  d4  SiJ  *  20  aeu 


(7) 


where  K  is  the  elastic  bulk  modulus  and  G  is  the  elastic  shear  modulus. 
These  moduli  can  be  functions  of  the  invariants  of  either  the  stress  or  the 
strain  tensor.  Accordingly,  it  is  assumed  that  K  -  K(Jt  ,  JJ  ,  JJ)  and 
G  -  G(J,  ,  JJ  ,  JJ)  or  K  -  K(I,  ,  II  ,  IJ)  and  G  -  G(I,  ,  H  ,  IJ)  . 
Equation  7  can  be  written  in  terms  of  the  hydrostatic  and  deviatoric  compo¬ 
nents  of  the  strain  and  stress  increment  tensors;  i.e., 


.  E 
d£kk 

1 

dJ 

3K(j, ,  JJ , 

JJ) 

< 

1 

dS 

2GU,,  JJ, 

JJ) 

(8) 


or 


dJ,  =  3K( I, ,  II,  IJ)  de^k 

dS^  -  2G(I,,  IJ,  IJ)  de^j  (9) 


where 


JJ 


de 

de 


E 

kk 

E 

ij 


11  “  ekk 
IJ  "  \  eij  eij 
IJ  ‘  3  ®ij  ®jk  eki 


increment  of  elastic  volumetric  strain 
elastic  strain  deviation  increment  tensor 
third  invariant  of  the  stress  deviation  tensor 
first  invariant  of  the  strain  tensor 
second  invariant  of  the  strain  deviation  tensor 
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third  invariant  of  the  strain  deviation  tensor 


In  order  not  to  generate  energy  or  hysteresis  within  the  elastic  range,  the 
elastic  behavior  of  the  model  must  be  path  independent.  The  material  should 
then  possess  a  positive  definite  elastic  internal  energy  function  W  which  is 
independent  of  stress  path.  The  strain  energy  function  can  be  written  as 


f  J 

W  -  J  0,  ,  de 


ij  ij 


r1J  .  1  .  .  s  T  ,  <*!■,  1 

l  (Sij  +  3  Jt  6ii)  9K(J,,  JJ,  JJ)  6ij  +  2007.  JJ,  JJH 


r‘  '  J*  JT  r1J  i 
”  J  9k( j,  ,  jj,  jj)  dJl  +  J  2007711717)  sij  dSij 


d(J, ) 


l  l8K(Jlf  JJ.  JJ)  +  l 


2g(j,,  jj,  jj: 


In  order  for  W  to  be  independent  of  stress  path,  the  integrals  in  Equation 
10  have  to  depend  only  on  the  current  values  of  J,  and  JJ  .  Therefore,  the 
bulk  and  shear  moduli  have  to  be  expressed  as 


K  -  K(  J, ) 


G  =  G( JJ) 


Similarly,  the  complementary  strain  energy  function  W  can  be  written  as 


w°  -  J 


'lj  doiJ 
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(12) 
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Substitution  of  Equation  1  into  Equation  12  leads  to 


’ij  e 


D  r  J  kk  P  P 

W  -  J  Hp  51j+eiJ)EK(I1,IMJ)(dekk-de^k)«  +  2GU,  ,IJ,IJ)(de  -de[  )] 


ij  ij,J 
(13) 


P  P 

But  dekk  and  de^  equal  to  zero  under  elastic  deformation.  Hence 
Equation  13  becomes 


D  r  ^  ekk 

Wu  -  J  (^«1J+eiJ)[K(Il,Ii,Ii)dekk6lj+2G(I1,IMJ)deij] 


2  IJ 

I  ~K  (I,,  IJ,  IJ)  d  [ek  ]  ♦  J  20  (I,,  IJ,  IJ,) 


Equation  14  indicates  that  for  vr  to  be  independent  of  the  strain  path,  the 
integrals  in  Equation  14  have  to  depend  only  on  the  current  values  of 
I,  and  IJ  .  Therefore,  the  bulk  and  shear  moduli  can  also  be  expressed  as 

K  -  K(ekk)  -  K(l, ) 

(15) 

G  -  G(IJ) 

It  can  be  concluded  from  Equations  11  and  15  that  the  bulk  modulus  should  be 
expressed  as  a  function  of  either  I,  or  J,  and  the  shear  modulus  should  be 
related  to  either  IJ  or  JJ  .  Further,  K  and  G  must  always  be 
positive.  Since  during  elastic  deformation  the  hardening  parameter  tc  is 
constant,  the  bulk  and  shear  moduli  can  also  be  expressed  as 

K  -  K( J, ,  k)  =  K(I, ,  k) 

(16) 

G  -  G(JJ,  k)  -  G(  I J ,  <) 
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Plastic  Strain  Increment  Tensor 

The  plastic  strain  increment  tensor  is  given  by  Equation  6  where  the 

loading  function  f  is  given  by  Equation  H  or  5.  The  hardening  parameter  in 

Equation  4  or  5  could  be  taken  as  being  equal  to  the  plastic  volumetric  strain 
P 


'kk 


thus 


'kk 


(17) 
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The  use  of  Equation  17  will  allow  the  cap  to  expand  as  well  as  to  contract 
(Figure  2).  The  plastic  loading  criteria  for  the  function  f  are  given  as 


K 


_3f 

9a 


do 


ij 


U 


>  0  for  loading 
-  0  for  neutral  loading 
<  0  for  unloading 


(18) 
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Because  de^  -  0  during  unloading  or  neutral  loading,  as  well  as  for 
f  <  0  ,  Equations  7  and/or  8  and  9  are  used  to  determine  the  purely  elastic 
strain  changes.  The  prescription  that  neutral  loading  produces  no  plastic 
strain  is  called  the  continuity  condition.  Its  satisfaction  leads  to  coinci¬ 
dence  of  the  elastic  and  plastic  constitutive  laws  during  neutral  loading. 

Like  the  elastic  behavior,  the  plastic  stress-strain  relation  can  be 
expressed  in  terms  of  the  hydrostatic  and  deviatoric  components  of  strain. 
Applying  the  chain  '-ule  of  differentiation  to  the  right-hand  side  of 
Equation  6  yields 
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Multiplying  both  sides  of  Equation  19  by  6^  gives 


(19) 
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The  devlatoric  component  of  the  plastic  strain  increment  tensor  de^  can  be 
written  as 


(21) 


Substitution  of  Equations  19  and  20  into  Equation  21  yields 


dl  9f  „ 

FV 7T  rVTT  ij 


(22) 


In  order  to  use  Equation  19  or  Equations  20  and  22,  the  proportionality 
factor  dA  must  be  determined.  This  can  be  accomplished  in  the  following 
manner.  From  Equations  5  and  17  the  total 

3f  hi  1  3f  p 

df  =  3^  +  ij 

I 

In  view  of  Equations  8  or  9  and  20,  Equation  23  becomes 


derivative  of  f  becomes 

asu  *  rr  d4  • 0  <23> 
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Substituting  Equation  1  into  Equation  24  results  in 
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D  P 

Substituting  the  values  of  dek(<  and  de^  from  Equations  20  and  22,  respec¬ 
tively,  into  Equation  26  yields 


3K  dekk  fjr +  4*  sij  deij " 9K  dX  (fS:)  +  G  dx  {rht) 
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Solving  for  dX  gives 


(28) 


Total  Strain  Increment  Tensor 

The  total  strain  increment  tensor  can  be  obtained  by  combining  Equations 
1,  7,  19,  and  28: 
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Similarly,  the  stress  increment  tensor  can  be  written  as 
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Equation  29  or  Equation  30  is  the  general  constitutive  equation  for  an 
elastic-plastic  isotropic  material.  To  use  either  of  these  equations  it  is 
necessary  only  to  specify  the  functic  .al  forms  of  K  ,  G  ,  and  f  and,  of 
course,  to  determine  experimentally  the  numerical  values  of  the  coefficients 
in  these  functions. 

Behavior  in  Tension 

The  tension  cutoff  is  triggered  in  this  model  whenever  the  following 
relation  is  satisfied: 


y  -  T  2  0  (31) 

J 1 

in  which  T  is  the  maximum  value  that  the  mean  hydrostatic  tension  y  can 
attain.  The  volumetric  strain  which  occurs  during  tension  is  computed  from 
the  following  relation: 
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EXAMPLE  MODEL  FIT  FOR  ISST  LAYER  2 

To  illustrate  the  new  model's  capability  for  fitting  iSST-type  material 
behavior,  selected  mathematical  forms  of  the  various  response  functions  were 
used  to  fit  the  material  properties  for  ISST  Layer  2  shown  in  Figures  ^ 
through  10.  In  the  following  model  response  functions,  tension  is  considered 
positive. 

Yield  Conditions 

The  yield  conditions  consist  of  an  ultimate  failure  surface  of  the  form 


f  -  F(J, ,  YJ J)  -  VjJ  -  [A  -  Cexp( BJ,  -  B1J,  )J 


and  a  work-hardening  cap  of  the  form 


f  -  H(Jj  ,  VJT  .  e£k)  -  (J,  -  L)2  +  R2Ji  -  (X  -  L)2  =  0 


(33) 


(3M) 


where 


Jt  -  first  Invariant  of  the  stress  tensor 


J2  ■  second  invariant  of  the  stress  deviation  tensor 
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T  (Figure  3). 
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A,  B,  B1 ,  and  c  »  material  parameters  (constants) 
p 

ekk  *  plastic  volumetric  strain 

P  P 

L(ek(<)  and  X(ekk)  "  respectively,  the  values  of  J1  at  the  center  of  the 

cap  (Equation  3*0  and  the  intersection  of  the  cap  with 
the  J,  axis,  and  are  related  as: 


l  -  X  -  R[A  -  Cexp(B£  -  Bit  )] 


(35) 


where  R  is  the  ratio  of  the  major  to  the  minor  axis  of  the  cap  and  is  given 
by: 
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where  RO,  R1 ,  R 2,  R3»  R4,  and  R5  are  material  parameters  (constants),  and 

(4(ekk^  if  1(ekk^  *  0 

L(ekkJ  “  \  (37) 

(°  lf4(ekk)>0 

Hardening  Function 

The  hardening  function  is  chosen  to  be 

ekk  *  W[exp(D.LL)-l~D.LL.exp(D1  .LL)]  -  D2.LL*.exp(D3. LL)  (38) 

in  which 


LL  -  L  - 


(39) 


where 

W,  D,  Dt ,  D2,  and  D3  =*  material  constants. 

W  in  Equation  38  defines  the  maximum  plastic  volumetric  compaction  that  the 
material  can  experience  under  hydrostatic  loading.  The  numerical  value  of 
ftft  is  the  solution  of  the  following  equation 

_  _  _ 2 

ftft  -  R[A  -  C.exp(B. ftft  -  Bl.lft  )]  -  0  ( 40) 


Bulk  Modulus 

The  bulk  modulus  was  taken  to  be  a  function  of  the  plastic  volumetric 
strain  and  the  first  invariant  of  the  stress  tensor 
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in  which 
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where  KI,  K1 ,  K2,  K3.  K4,  K5,  K6,  K7t  K8,  and  RKJ  are  material  parameters 
(constants) . 

Shear  Modulus 

The  shear  modulus  was  taken  to  be  a  function  of  the  plastic  volumetric 
strain  and  the  second  invariant  of  the  strain  deviation  tensor: 
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parameters  (constants).  SL  and  eekk  in  equation  47  are  given  by 
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where  G1 3  i3  a  material  parameter  and  the  initial  value  of  SL  is  zero. 
Numerical  Values  of  Fitting  Parameters 

The  above  model  includes  a  total  of  39  fitting  parameters  (or  material 
constants);  the  numerical  values  of  these  parameters  used  to  fit  the  material 
properties  specified  for  ISST  Layer  2  are  given  in  Table  1.  Figures  11-21 
depict  the  stress-strain  and  strength  behavior  predicted  by  the  model;  the 
plot  scales  permit  direct  comparison  with  the  material  property  plots  given  in 
Figures  M-10. 
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Table  1.  Numerical  values  of  fitting  parameters  for  ISST  Layer  2. 


Model 

Fitting  Parameters 

Numerical 

Name 

Notation f  Unit 

Value 

Bulk  Density 

p  ,  g/cm^ 

1.90 

A  .  MPa 

173.2051 

failure 

B  ,  (MPa)"' 

0.0013 

Envelope 

_7 

Parameters 

B!  ,  (MPa)  ‘ 

12  x  10  ' 

C  ,  MPa 

172.916JI 

RO  ,  -- 

2.5 

R1  .  — 

0.0 

R2  ,  (MPa)-’ 

0.0 

R3  .  — 

2000.0 

Hardening 

R9  .  — 

1.5  X  10* 

field 

R5  ,  (MPa)-1 

0.02 

Surface 

Parameters 

W  ,  — 

0.22 

D  ,  (MPa)-1 

0.0105 

> 

Dl  ,  (MPa)-1 

0.095 

D2  ,  (MPa)-2 

0.00015 

D3  ,  (MPa)-' 

0.09 

K1  ,  MPa 

2000.0 

K1  .  “ 

0.9 

K2  .  — 

2.0 

K3  ,  MPa 

105 

Bulk 

K8  ,  - 

2000.0 

Modulus 

K5  ,  — 

29.3 

Parameters 

K6  ,  — 

0.8 

K7.  ,  “ 

105 

K8  ,  — 

15.0 

RKJ  ,  MPa 

-1500.00 

GI  ,  MPa 

1500.0 

GI  ,  — 

0.9 

G2  ,  — 

2.0 

G3  .  MPa 

0.75  x  105 

G9  ,  -- 

2000.0 

G5  ,  -- 

22.5 

Shear 

G6  ,  — 

220.0 

Modulus 

Parameters 

G7  ,  — 

900.0 

G8  ,  — 

0.219829 

G9  ,  — 

0.063 

G10  ,  — 

12000.0 

Gil  . 

0.163 

GI  2  ,  — 

120.0 

Gl 3  .  “ 

0.185 
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Figure  1.  Uniaxial  strain  (UX)  stress  paths  for  ISST  materials 
(extract  from  Reference  2). 


MSS 

rcMXd 


■*'s. 

'.‘iL'i! 

J^TST*5 

%r«BTV* 
IV  v •• 

cy>& 


£&&& 

%;,i 

vJv 

■*  \  •  _■ 


■V/v 

Cv'vs'a 

Ir^ 


wv* 

*  V  rf"  V 

$1 

•Crab 


VERTICAL  STRAIN,  PERCENT 

Figure  5.  Material  properties  for  ISST  Layer  2:  uniaxial  strain 
compressibility  relation  to  =  800  MPa. 
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Figure  10.  Material  properties  for  ISST  Layer  2:  triaxial 

compression  stress-strain  relations  for  high  confining 
pressures. 
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Figure  13.  Strain-dependent  cap  model  fit  for  ISST  Layer  2 
uniaxial  strain  compressibility  relation  to 
o  ■  860  MPa. 
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Figure  17.  Strain-dependent  cap  model  fit  for  ISST  Layer  2: 

UX  stress  path  to  o  ■  2,500  MPa  and  TX  failure 
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Figure  18.  Strain-dependent  cap  model  fit  for  ISST  Layer  2: 

triaxial  compression  stress-strain  relations  for 
low  confining  pressures. 
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Figure  20.  Strain-dependent  cap  model  fit  for  1SST  Layer  2 
triaxial  compression  principal  stress 
difference  -  vertical  strain  relations  for  low 
confining  pressures. 
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